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AGN. II: X-ray excited funnel flow 

A. Dorodnitsyn^'^, T. Kallman^, and D. Proga^ 
ABSTRACT 

We have calculated a series of models of outflows from the obscuring torus 
in active galactic nuclei (AGN). Our modeling assumes that the inner face of 
a rotationally supported torus is illuminated and heated by the intense X-rays 
from the inner accretion disk and black hole. As a result of such heating a strong 
biconical outflow is observed in our simulations. We calculate 3-dimensional 
hydrodynamical models, assuming axial symmetry, and including the effects of 
X-ray heating, ionization, and radiation pressure. We discuss the behavior of a 
large family of these models, their velocity fields, mass fluxes and temperature, 
as functions of the torus properties and X-ray flux. Synthetic warm absorber 
spectra are calculated, assuming pure absorption, for sample models at various 
inclination angles and observing times. We show that these models have mass 
fluxes and flow speeds which are comparable to those which have been inferred 
from observations of Seyfert 1 warm absorbers, and that they can produce rich 
absorption line spectra. 

Subject headings: acceleration of particles - galaxies: active - hydrodynamics 
-methods: numerical - quasars: absorption lines - X-rays: galaxies 

1. Introduction 

One of the insights provided by observations of Seyfert galaxies and some quasars is 
the prevalence in their X-ray spectra of spectral lines and bound-free continua from ions 
of intermediate-Z elements . Early observations of Seyfert 1 galaxies using proportional 
counters and solid state detectors revealed spectra with strong absorption features in the 
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0.1-10 keV range (iHalperru Il984j ). These features were attributed mostly to the edges of 
hydrogen and hehum - hke oxygen. The term "warm absorber" was proposed owing to 
the fact that the observed X-ray absorbing gas has an electron temperature lower than it 
would be if a similar level of ionization were produced by coUisional ionization. However, 
more detailed spectroscopic studies were hampered by the limited X-ray resolution of the 
ASCA and ROSAT satellites. The grating spectrographs on the X-ray telescopes Chandra 
and XMM-Newton provide unprecedented spectral resolution up to ~ 10 keV. These show 
that X-ray spectra obtained from ~ half of low-red-shift active galactic nuclei (AGN) contain 



(Kaspi et al. 


* — 
2002; 


SteenbruEcee 


2005) 



absorbing g as has been con f irmed in the majority o f AGNs which are bright enough to allow 
detections ( Reynold^ 1997 : McKernan et al. 2007 ) . The re is also a partial correspondence 



between UV and X-ray absorbers (jCrenshaw et al.lll999l ). 



X-ray observations of warm absorbers are consistent with the Seyfert 1/Seyfert 2 di- 
chotomy. For example, the properties of the X-ray emission in the Seyfert 2 galaxy NGC 1068 
corresponds to the scatte red emission expected from warm absorbers in Seyfert 1 galaxies 
jKinkhabwala et all 12002 ). 



Constraints on the position and dynamics of the X-ray absorbing gas can be deduced 
from the observed widths and virial argument s, and also from the variability studies of these 
spectra (IBehar et al.ll2003l : iNetzer et al.ll2003l ). These show an absence of correlated response 
of the warm absorber gas to rapid changes (~ days) of the continuum. This implies that the 
ionization time scale in the warm absorber gas is long (> months). Combined together, the 
line blueshifts, widths, and time variability analysis favors an origin of the warm absorber gas 
at -R > 1 pc away from the BH. This estimate coincides wit h the likely location of abs orbing 
matter responsible for obscuration in Seyfert 2 galaxies (IKrolik fc BegelmanI Il988 ). The 



existence of an outflow from the torus has been s uggested by 



Krolik fc BegelmanI (11986 



1988), and as the source of warm absorber flows by Krolik fc KrissI ( 19951 . 2001). 



It is believed that this matter is in the form of a molecular torus which is responsible for 
obscuring the broad line region in Seyfe rt 2 galaxies, and which is thought to exist in most 
low and intermediate luminosity AGN (lAntonucci &: Millerl Il986l ). A growing body of di- 
rect observational evidences advocates for the existence of the obscuring torus. Mid-infrared 
high spatial resolution studies of the nucleus of NGC 1068 using the Very Large Telescope 
Interferometer ha ve resolved a dusty structure which is 2.1 pc thick and 3.4 pc in diameter 
(jjaffe et al.l 12004 ). Observations support a multi-temperature model: the temperature of 
the warm component was established to be 300 K and inside of it a second, compact and hot 
(>800K) component has been found. Further studies of NGC 1068 systematically reduced es- 
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timates of the temperatures of different components (jPoncelet &: Perrinll2006l ). Observations 
of the Circinus galaxy, which is among the closest prototype Seyfert 2 galaxies, also revealed 
a dense and warm T > 300 K conaponen t at about 0.2 pc from BH and cooler T < 300 K 
component at Ipc (ITristram et al.l 120071 ). If the hotter component is located closer to the 
X-ray source, it may be attributed to the inner part of the torus, heated by the radiation 
of the compact nucleus. Although the evidence is strongest for nearby active galaxies, there 
is also a strong motivation to think that within the same obscuring torus paradigm exist 
those quasars whose central regions are heavily obscured by gas an d dust (Type II quasar s). 
Evidence for this comes from spectro-polarimetric observations by IZakamska et al.l (120061 ). 



This paper is part of a series whose main goal is to test the hypothesis that the torus is 
the origin for the warm abso rber flow. Preliminary results of this work have been reported 



Dorodnitsyn et al.l (120081 ) (Paper 1), in which we presented the results from a sample 



m 

model and showed that the adopted model is promising in explaining the warm absorber 
phenomenon. In this paper we provide more details of our methods, and display results of 
models which span the space of input parameters. We present and discuss the hydrodynamic 
quantities which characterize our models: mass fluxes, velocity fields, and temperature struc- 
ture. We also show sample X-ray spectra, which we will discuss extensively in a later paper 
of this series. 

Our approach can be described as having three basic parts: i) setting up initial condi- 
tions, which requires defining an initial torus configuration and making assumptions about 
the external source of radiation; ii) implementation of the wind driving force (local heating- 
cooling rates and radiation pressure force) and actual 2D hydrodynamical calculations. The 
latter includes the numerical solution of the time-dependent 2D (so called 2.5D) system of 
equations, which takes into account centrifugal forces, and radiation pressure and heating 
terms; iii) calculating of the X-ray line spectra using a code which adopts Sobolev radiation 
transfer and ionization calculations for plasma in the intense X-ray field. Each of these steps 
is described in what follows. 



2. Governing equations 

We solve the following system of equations: 



| + V.(pv) = 0. 
p(| + (v.V)v 



-Vp-pV$ + pgrad, 



(1) 

(2) 
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- + V-(v(e + p + p«l>)) = i/. (3) 

These are the conservation equations for: mass, momentum and energy. Heating and coohng 
processes are described by the function iJ(ergcm~^ s~^); e - is the sum of the kinetic and 
internal energy densities: e = + e. These equations should be supplemented by 

the equation of state which we assume to be polytropic: P = K p'^, where 7 = 1 + 1/n, 
and n is the polytrope index and P = (7 — 1) e. A one-component, one-temperature T = 
Pp/ pTZ, where p is the mean molecular weight per particle, TZ = 8.31 ■ lO^erg g~^ is the 
universal gas constant and plasma with 7 = 5/3 is assumed to constitute the flow. All three 
components of the flow velocity v = {vr, vg, v^) are calculated, assuming azimuthal {d/dcp = 
0) symmetry. Equations ([I])-© are cast in a non - dimensional form with the characteristic 
scales set by the properties of the plasma orbiting at a characteristic distance, Rq from 
a black hole of mass Mg ( in units of IO^Mq). The characteristic scales are respectively: 
to = -Rq^^/ V GM ~ 4.5 ■ 10^^ Tpc^ Mq^^'^ (s) for the time, where r^c is the distance in parsecs, 
and Vo = ^GM/Rq ~ 6.6 ■ WmI'"^ r^c'"^ (cms'^) for the velocity. 



3. Forces driving the flow 

Heating and cooling of the gas 

The forces accelerating the wind in our model result from the gradient of gas pressure 
and from radiation pressure. The thermodynamic properties of X-ray heated gas depend on 
the spectrum of the incident radiation as well as on the local atomic physics. Under the 
assumption of photo-ionization equilibrium the thermodynamic state of photoionized gas 
can be pararneteri zed in terms of the ratio of radiation energy density to baryon density 



(ITarter et al.lll969l ): 



e = 47rF,/n, (4) 
where Fx = L^e~'^ / (47rr^) is the local X-ray flux, Lx is the X-ray luminosity of the nucleus, 

and T = Kpdr - is the optical depth, and n is the number density. We assume that the 
Jo 

attenuation is dominated by Thomson scattering n = 0.2(1 + Xh) — 0.4cm^g~^, where 
Xn is the mass fraction of hydrogen, and the factor e~^, accounts approximately for the 
attenuation of the radiation flux on the way from the source toward a fiducial point. The 
methods adopted in t his paper fo r treat in g the e f fects of radiation are essentially the same 
as those described in iProga et all (I2OOOI ): iProgal (120071 ). and have been applied to various 
problems in the study of AGN and X-ray binaries. 
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Assuming that the there is a fraction /x of the total accretion luminosity Lbh available 
in X-rays and that the disk radiates a fraction F of its Eddington luminosity Lgdd = 1.25 ■ 
10^^ Mg we estimate: ^ ~ 4 ■ 10^ ■ /x F M6/(A'"23 rpc), where A^23 is the column density in 
10^^ cm~^. If the dynamical time within the flow is much larger than the characteristic 
time of the photoionization and recombination then the ionization balance is determined by 
the condition of photo-ionization equilibrium. The rates of Compton and photo-ionization 
heating and Compton, radiative recombination, brems strahlun g and line cooling are then 
given by approximate formulas, modified from those of iBlondinI (119941 ). for these processes: 



Fic(ergcm ^s 
for the Compton heating - cooling; 



.9- 10"''''^(Tx-4T), 



(5) 



Fx(ergcm-3s-i) = 1.5 • lO^^i ^i/4 5^-1/2 _ ^^^-i^ 

for the photo-ionization heating-recombination cooling , and for the bremsstrahlung and 
line cooling: 



A(ergcm-^s~^) = 3.3 ■ lO^^^T^/^ 

+ (4.6 ■ 10-1^ exp(-1.3 ■ io5/T)e(~o-«-°-9^")T-l/2 + lO^^^) 5. (7) 



These formulae have been originally derived for a 10 keV bremsstrahlung spectrum (Tx = 
2.6 ■ 10^ K) and were fo und to be in a reasonable ( ~ 25%) agreement with numer ical sim- 
ulatio ns (IBlondinI I1994J ). Equations ([3])-([7j) are slightly modified versio n of those of iBlondin 
( 19941 ) ■ which accommodates new atomic data. Using the XSTAR code ( iKallman &: Bautista 
2OOII ) we recalculated heating-cooling rates for the incident spectrum which is a power law 
with energy index a, and found results essentially equivalent to those given by equations 
©-([Tj). Notice that in the case of a bremsstrahlung spectrum a formal value of a = should 
be used in ([7]). For a power law with energy index a = 1.1 the results differ by < 30 % (see 
Figure [1]). Given these rates of energy deposition from the radiation to the flow, we write the 
total radiative heating-cooling function: H = Tic -|- Fx — A. We have also performed several 
runs of our hydrodynamical models with different assumptions about heating-cooling, and 
found no ir nportant differe nce in the flow dynamics if using equations ([3])-([7j) or the original 
formulae of IBlondinI (119941 ). and also between bremsstrahlung and power law spectra for sev- 
eral values of a. It appears that, for example, the effects of the optical depth are much more 
important. That is, the difference between curves for the power law and the bremsstrahlung 
spectrum at small ^ (correspondingly high density) in Figure [1] becomes unimportant. 
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Fig. 1. — Comparison of the results from the XSTAR X-ray photo-ionization code and 
analytic approximations ([3])-([7j). Curves are plotted for different values of ^ ranging from 
^ = 1 (lower curves) to ^ = 10^ (upper curves). Vertical axis: total radiative heating-cooling 
function H = Tic + Fx — A in units of 10~^^ ergcm^s"^. Horizontal axis: log(temperature). 
Curves; solid: XSTAR; dot-dashed: analytic. 
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The radiation pressure 

The radiation pressure force consists of the force due to continuum absorption gfcont 
F\jyk/c and due to hues: 



^rad = (FuvK/c)M(t), 



where M{t) is the force multiph er (ICastor et al.lll975r) . and Fyv is the local UV flux. We 
make use of the particular form (jOwocki et al.lll988l ): 



v(l-a) 



1) /r^-" 

J / max ' 



(9) 



M{t) = kt-'^ ((1 + r, 

where t = r/rj is the optical depth parameter, rj = ni/cTe is the line strength parameter, 
(Te is the Thomson cro ss-se ction, and Tmax = trima^. A parameter T^max was introduced 
by lOwocki et al.l (Il988l ) and IStevens &: KallmanI (jl990l ) in order to limit the effect of very 
strong lines. That is, they assume a line number distribution which satisfies: dN/ [drj dv) ~ 
7]""^ exp(— 77/?7max), where N{rj, u) is the line number dist ribution. If r? r nax — » oo, so that lines 
are distributed as a power law, one recovers the result of ICastor et al.l (119751 ): M{t) ~ kt~"'. 
In the opposite case of Xmax 0, the force multiplier is independent of t, and Mmax ~ krj^^^. 
As a result of this maximum line strength cutoff a correction factor appears in the relation 
for M(t), ([9]). The dependence of k and 77max on £ has been numerically calculated and then 
fitted by the analytical formulae ( Stevens fc KallmanI ll99ol ): 

= 0.03 + 0.385 exp(-1.4^°-6), (10) 



logio Vn 



6.9exp(0.16^°-^), logio^<0.5, 
9.1 exp(-7.96 ■ lO^^ ^), log^Q ^ > 0.5. 

From these, one can see that M{t) can depend sensitively on the ionization parameter. 
Taking a fiducial a = 0.5 (the value, adopted in all our calculations), one finds that Mmax = 
585 at ,^ = 0, then has two local maxima: Mmax = 724 at ,^ = 0.3, and Mmax = 743 at 
^ = 3.1. Mmax then drops to 1.7 at ,^ = 100 and decreases gradually to Mmax = 0.01 at 
^ = 1000. 



3.1. Initial configuration: rotating torus with arbitrary Compton optical 

depth. 



We begin from a rotating toroidal configuration which is in equilibrium in the external 
gravitational field of the BH. The equation of state of the torus interior is described by the 
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polytrope P = Kp^^^l"^ . The distr i bution of the density (or pressure) in the torus interior 



was given by iPapaloizou fc Pringld (119841 ) (in what follows refer to PP-torus for short), who 
assumed that the distribution of the specific angular momentum inside the torus is constant. 
In our case such a torus would not be in equilibrium becaus e of the radiation pres s ure fr om 



the central object. Thus we modify equilibrium equations of IPapaloizou fc Pringld (119841 ) to 
include the radiation pressure term. Since this cannot be done in a closed analytical form, 
we can write an approximate equation: 



V 



1 _ r p-^{'^) 

C\. (11) 



p n + l\ r 2r2 sin^ Q 

Note that equation ([TT]) must be understood as a bridging formula between two limiting 
cases: optically thin e""^ ~ 1 (in which case it is the PP-torus with 1 — P reduced gravity), 
and optically thick, when (yfrad ~ (PP-torus case). A constant C in the equation (fTTj) . 
parameterizes the distribution of the torus models and is connected with the distortion of 
the torus (this is described in more detail below). Including the radiation pressure reduces 
the effective gravity, and thus the torus gas needs less entropy to sustain it against vertical 
collapse. In both of these limiting cases this equation is exact. 

Notice that the problem of toroidal equilibrium in the presence of heating (or other 
radiation transfer effects) introduces a characteristic length scale through the optical depth 
r, leading to non-self-similarity of the model. Equation (fTTj) was derived by assuming that 
the distribution of the specific angular momentum inside the torus is constant. Choosing 
non-dimensional units and working in terms of zu, the cylindrical radius in units of i?o, if 
we define the non-dimensional density p such that piw = 1) = 1, and the non-dimensional 
pressure P and internal energy e such that: P = (7 — l)e, and e{w = 1) = Cq then 

Ti (1 r) (\ 1 \ 

en = -, I 7: — C\ . The inner and outer edges of the torus are lo- 

cated at w~ and cc"*", respectively. Bounded configurations exist only for < C < 0.5 and 
the distortion of the torus is described by the parameter d = {w~ + ti7+)/2 = 1/ (2C). The 
boundary of the torus is matched to the exterio r by the condition P = 10~ ^. The PP-torus 



is unstable to non-axzsymmethc perturbations (IPapaloizou fc Pringld 1 19841 ). However, this 
effect cannot be numerically investigated in the azimuthal symmetry which we adopt, since 
no signals can propagate in the direction. At w > 1, matter that constitutes the torus 
has an excess of angular momentum with respect to the local 'equilibrium' keplerian value, 
z) > l{r), and vice versa, in the inner parts of the torus w < 1, l{w, z) < l{r). It is the 
internal pressure of the torus (equation ITTl) that inhibits matter from settling to smaller (or 
larger, depending on angular momentum) orbits. The gas first evaporates from the part of 
the torus which is closer to the source of radiation and tends to settle at larger zu, as soon 
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as the back pressure supporting it drops. 



We begin our simulations from the stationary configuration determined from equation 
f ITT]) . We follow the torus evolution as it is being heated by X-rays. No replenishing of 
the gas which constitutes the initial torus is provided: Therefore, the torus will eventually 
lose all its mass and will completely evaporate. However, in the regime we are looking 
for, the evaporation is not dramatic and does not significantly deplete the torus during the 
characteristic dynamical time. 

In the following sections we will show that the existence and character of the fiow from 
the heated torus depends critically on the geometry. That is, it depends on the divergence of 
the fiow streamlines, the strength and incident angle of the X-ray illumination, and on the 
direction of the effective gravity in the rotating frame of the torus. The flow is intrinsically 
two-dimensional, and t herefore cannot be ad e quate ly described a priori by ID models, such 
as those preformed by IChelouche fc Netzerl (120051 ). Furthermore, the shape of the torus, 
and thus the launching surface for the flow, is affected by the flow. So the tor us interior 
cannot be considered as a boundary condition (e.g. as in 
to include it in the computational domain. 



Balsara fc Krolik ( 19921 )): we need 



4. Methods 

For our computations we adopt a spherical-polar coordinate system {r,9), extending 
the computational domain {ri,9j} from = 0.01 to rout = 50 in radius, and from to 
TT in the polar domain making no assumption about equatorial symmetry. The number of 
points in the radial, N^, and polar, Ng, directions are taken to be equal: = 140, in low 
resolution and = 300 in high resolution grids. The {rj} grid is non-uniformly spaced, i.e. 
r2 = r, + (rout - n^){kl^^'''-'^ - l)/((fc^/(^^-^) - 1)), and r,+i = r, + (r, - r^.^^y^^^"^), 
for i = 2, Nr — 1, and the reflnement factor is kr=4:. In order to achieve better resolution 
of the flow itself rather than the torus interior we also adopt a polar grid which has non- 
uniform spacing 56i = 6i — Oi_i, so that the maximum reflnement is approached at 6 = tt/4: 

= 5i/k]'^^'~^^ at < ^ < 7r/4, and 6 = n/A: 5i+i = 5ik]'^^'~^'^ at Ti/A<e < tx/2 (and 
analogously spaced in the southern hemisphere) . Boundary conditions are axially symmetric 
at 6* = 0, vr and outflowing at rin and rout 

To so lve numerically the sys tem of hydrodynamical equations ([I])-([2]) we use the code 
ZEUS2D Jstone fc Norma Jl992h . Note that the characteristic time of X-ray heating/ cooling 



can be much shorter than the dynamical time, which in such a case introduces strong stiff- 
ness to the system of equations ([I])- ([3]). To overcome this difficulty, some modiflcations have 
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been made to the code. The most important one is the implementation of a fully implicit 
update of the energy in equation ([3]) just prior to the transport step in ZEUS2D. Addionally, 
we account for the radiation pressure (equation ([S])) term. As an initial test we have evolved 
a toroidal distribution of matter for two rotational periods and found the configuration to 
be stable. The gas is illuminated by the incident X-ray radiation with a power law spectrum 
with an energy index, a = 1. The heating/cooling rates are described by the approximate 
analytical formulae give in equations ([3])-(j7]). 



Warm absorbers 



We test output of our hydrodynamical models against the ability to predict warm ab- 
sorber spectra. To do this we use the output from the hydrodynamical code, of p, v, and T as 
an input to the calculation of X-ray line and photoelectric absorption spectra. The numeri- 
cal code has been specifically developed for calculati on of spectra in the X-ray domain and 
makes use of procedures developed for the XSTAR (iKallman fc Bautistall200ll ) code, while 
calculating the ionization structure an d distribution of opacities, and treating the radiation 
transfer in the Sobolev approximation (IRybicki fc Hummerlll983l ). Although the goal of this 
paper is to show that pure hydrodynamic 2D models can produce warm absorber spectra, we 
present here only sample spectra, assuming pure absorption. We postpone a more detailed 
discussion, including a full 3D transfer calculations, to a separate pubhcation. 



5. Results. 

The most important parameters which determine the properties of the warm absorber 
flow are the initial Compton optical depth = t{6 = 90°) of the torus (or equivalently the 
maximum initial torus density rimax)? and the distance from the BH, Rq. We also explore 
the dependence on F and d. Other parameters are chosen having some typical values: the 
mass of the black hole: Mbh = 10^ Mq, the Compton temperature of the X-ray radiation 
Tx = 10 keV, and the fraction of X-rays and UV radiatio n = fuy = 0.5. (For rotating 



flows exposed to a multi-temperature radiation, see e.g. iProga et al.l (l2008l )). The latter 



is con sistent with typical energy distributions of the radiation close to the BH (ILaor et al. 



19971 ). We neglect any changes in the BH luminosity. The important thermal time scales 



within the flow, namely the Compton heating and cooling time, tx and the dynamical time, 
tdyn, may be of the same order ~ tdyn ~ 10^° s. This is discussed in more detail later 
in this section. Thus the outflowing gas may not be in thermal equilibrium and adiabatic 
losses are likely to be important. Notice that a nearly hydrostatic Compton heated corona 
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can exist only at r < Ric = GMbr fi/TZT^ ^ 8- IO^^^Mq/T^j) cm, where T^j is the Compton 
temperature in terms of lO^K. In all of our models the major flow is located at r >> Rjc- 

We have calculated 20 models, including combinations for: = 1.3 (models Ai), and 40 
(models Bi); Rq = 0.5, 1, 1.5; F = 0.1, 0.3, 0.5 all with d = 2.5; and two models with d = 5 
(models Ci). These are summarized in Table 1 where some of the characteristic results from 
the computed models are presented. In what follows we describe in detail the cases which 
best illustrate the most important results. We also discuss the dependence of our results on 
parameters, based on the behavior of the ensemble of models. 

The model Aq is similar to that described in Paper 1, although the initial torus in the 
model which is described here has a different distribution of p and e (see equation (fTTj) . and 
the discussion thereafter), and smaller r^. In Paper 1 this model has been described in 
detail. Calculations presented here reveal more details and confirm the conclusions of Paper 
1. We begin here by describing results from model Bq, and later discuss how it differs from 
model Aq. 

Model, Bq has = 40, Rq = 1 and F = 0.5 and corresponds to a Compton thick 
(r^ ~ 40) torus having large rimax = lO^cm"^ and mass Mtor = 9.3 ■ 10^ M©. Results 
are displayed in Figure [2l where the evolution of the distribution of density is shown as a 
function of time (the density scale is such that corresponds to lO^cm"'^); Figure [3] where 
the distribution of pressure is shown at t = 3 (the pressure scale is such that corresponds to 
4.7- 10~^dyn/cm^) ; Figure Hlfor various quantities as a function of the inclination 6; Figure 
|5] (left panel), where the effect of the distortion parameter d is demonstrated; and Figure 
|6] showing horizontal 'slices' of the velocity and temperature at constant height, z. In the 
case of this model, the torus column is high enough to effectively screen the torus interior 
from penetrating X-rays. This leads to formation of a nearly pure funnel flow, i.e. the torus 
interior, and hence the shape of the surface responsible for launching and collimating the 
flow, is essentially unaffected by X-ray heating on time scales < trot- 

Here and in what follows we discuss the time evolution of our models in terms of t, 
measured in units of the characteristic time of rotation, to- After t = 1 a. high pressure 
region created by X-ray heating extends to r ~ 4.5 — 5 pc throughout the area that is not 
shadowed by the high density torus. At this time the torus is located at 6* < 50°. The 
distortion parameter has a value d ~ 2.5, i.e. the torus shape is almost unchanged from its 
initial value. This is shown in the upper left panel of Figure [2l Within the part of the flow 
which is not shadowed by the torus, high temperature gas expands in a spherical bubble with 
radius, r < 5.2 pc in which the temperature is T ~ 3 — 10 T^i^{r), where T^ir = 2.6-10^ M^/r^c 
is the local virial temperature. An axisymmetric region exists between < 0.75 pc and z <2 
pc where the temperature, T ^ lOTvir(r) . That is, high temperature, T ~ 3 ■ 10^ K, but 



- 12 - 



low density gas fills the torus funnel. The ionization parameter (equation Hj) in this region 
is ~ 10^ — 10^. The outer edge of the torus extends to ~ 4.25 pc in temperature, and to 
~ 4.5 pc in density contours. 

Figure [2] (lower left) shows density and velocity fields for model Bq at t = 3. Figure 
|3] shows that a high pressure region expands to height z ~ 6 pc from the equatorial plane. 
The torus inner edge is inferred from the temperature and density maps to be w~ ~ 0.83 pc. 
Inside the of this radius, which we refer to as the torus throat, the temperature is T ~ 
10^ — 10'^ K. A wide nozzle with (-CiJmax — '^min)/^max — 2.12, where 2;max — 0.4 pc is formed, 
having inner radius of ~ 0.85 pc. The torus outer edge is slightly shifted to zu~^ ~ 4.5pc. 
The values of ^min (the minimum ionization parameter along a radial line) and the column 
density vary significantly with the inclination angle. Figure H] shows the distribution of 
radial and poloidal velocity, ^, density and the rate of growth of number density with radius 
as function of ^ at t = 3 for model Bq. Near the axis, ^min(^ — 20°) = 10^ and the 
column density is N23 = 10~^. Note that if ^ tdyn, i-e. the gas is not in thermal 
equilibrium, then ^ is not as meaningful as when « tdyn- When tx > tdyn adiabatic losses 
strongly affect the temperature of the gas. At larger 6, the ionization parameter decreases: 
^jnin(6' — 25°) = 3-10^ and at higher inchnation, ^ gradually reduces from ^min(6' = 45°) = 12, 
eventually becoming ^min(6' — 60°) = 2.5. At a critical angle, 9 ~ 40°, a strong rise of the 
column density reflects the fact that the line of site penetrates the dense torus body rather 
that through the wind (c.f. Figure HI lower right). The column density increases, from 
A^23 = 0.3 at 6* = 45° to N23 ~ 100 at ^ 60°, providing total obscuration. Figure [3] 
also shows the position of the sonic surface determined by the relation Vp/cg = 1, where 
Vp = (f^ + VqY^'^ is the poloidal velocity and Cg = (7^T//i)^/^ is the speed of sound. Behind 
the torus a low entropy region exists which is bounded from the sides by a quasi- stationary 
shock. The existence of this structure can be understood from the following considerations. 
If the flow were perfectly symmetric in both hemispheres, then it should have = at 
z = 0, and the z = plane would be the equivalent of a rigid wall (reflecting boundary). 
Thus, if f z < behind the torus the formation of a shock structure is anticipated. Generally, 
this is the kind of picture one expects to observe from a supersonic wind flowing over a rigid 
obstacle. 

At t = 5 in model Bq (Figure [2], lower right), the density maximum is located at 
ti7 ~ 2 pc. The inner edge of the torus does not shift significantly from the position it has at 
t = 3: -uj^ ^ 0.75 pc in density maps (and ~ 1 pc in temperature maps); the outer edge is 
at ~ 4.3 pc. The temperature of the torus interior is in the range 10^ — 6 ■ lO'^K. A hot 
flow is located near the axis, bounded from the sides by the torus throat, and having high 
temperature: ~ few ■ 10^ K. A significant drop of ionization parameter ^ from ~ 6 ■ 10^ to 
~ 6, occurs again at 6* > 45° — 50° (c.f. Figure H]) , where the column density also rises from 
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N23 = 0.04, to N23 = 30 at 9 > 60°. The asp ect ratio of the torus is: A = Rn/H ry 1 i n 



accord with what is inferred from observations (iKrohk fc Begelmadll986l : iJaffe et aLll2004j ). 
At low inchnations, 9 < 10°, everywhere in the wind the poloidal component of the velocity 
is determined by Vr- However, at 6* > 50° inside the torus throat, the vg component is 
important, i.e. ve ~ Vr aX ru < Ipc. 

Model Aq has = 1.3, i?o = 1, T = 0.5 and is very similar to the model described 
in Paper 1. It differs from model in that the smaller optical depth of the torus interior 
cannot shield the gas from a extensive X-ray heating and the torus loses mass from large 
parts of its surface. The initial maximum density of the torus is n^aax = 10® cm~^ corresponds 
to initial torus mass, Mtor = 9 ■ lO'* Mq. Figure [7| shows the distributions of poloidal velocity, 
^, density and the rate of growth of number density with radius as a function of ^ at t = 3 for 
model Aq (in the same format as Figure H]). During the evolution, a region of high pressure 
extends from r ~ 4.5 — 5 pc at t = 1 to r ~ 12 pc at t = 3, and to r ^ 20pc at t = 5. The 
inner edge of the nozzle shifts slightly from uj^ ~ 0.8 pc at t=l to ru^ ~ 0.83 pc at t=3, and 
w~ ~ 0.75 pc at t=5. At later times the behavior of the model Aq is similar to models Ai 
and A^ , and can be inferred from Figure M 

It has been mentioned that in model Bq much of the torus interior is opaque to pene- 
trating X-rays. Remarkably, the minimum nozzle cross-section doesn't change much at late 
times, implying that the mass-loss rate becomes quasi-saturated. Note that in the case of a 
ID flow M is roughly set by the position of the sonic point, which in turn is set by gravity. 
In the case of a 2D nozzle, the mass-loss rate is determined by X-ray heating, gravity and the 
minimum nozzle cross-section. In the case of model Bq the latter remains almost unchanged 
in time. We believe this model is probably most representative in showing the key features 
of X-ray excited flow. However, models Ai may generally have broader angular patterns in 
which a warm absorber spectrum is observed, as will be discussed below. Only comparing 
synthetic spectra with observations can answer the question of what model is more adequate 
in describing the phenomenon of warm absorbers. 
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Model 




lid 


r 


d 


max,t=3 


^max,t=3 


^max,t=5 


^max,t=5 


M 


L — 


Mt-s 


Ai 


1.3 


0.5 


0.1 


2.5 


516 


155 


624 


332 


4.09 


10"^ 


6.54 ■ 10-3 


A2 


1.3 


0.5 


0.3 


2.5 


710 


317 


847 


330 


1.48 


10-3 


4.31 ■ 10-3 


A3 


1.3 


0.5 


0.5 


2.5 


707 


267 


760 


291 


2.34 


10-3 


2.14- 10-2 


A4 


1.3 


1 


0.1 


2.5 


547 


189 


514 


217 


1.76 


10-3 


1.66 ■ 10-2 


As 


1.3 


1 


0.3 


2.5 


526 


179 


605 


343 


9.65 


10-3 


6.34- 10-2 


Ae 


1.3 


1 


0.5 


2.5 


570 


235 


670 


337 


2.02 


10-2 


1.8 ■ 10-2 


A7 


1.3 


1.5 


0.1 


2.5 


360 


197 


413 


230 


6.21 


10-3 


1.64- 10-2 


As 


1.3 


1.5 


0.3 


2.5 


388 


169 


540 


310 


1.40 


10-2 


5.68 • 10-2 


A9 


1.3 


1.5 


0.5 


2.5 


317 


207 


663 


370 


2.66 


10-2 


1.23 • 10-1 


Bi 


40 


0.5 


0.1 


2.5 


673 


318 


522 


320 


3.56 


10-3 


7.38 • 10-3 


B2 


40 


0.5 


0.3 


2.5 


590 


257 


1004 


471 


1.16 


10-3 


1.87-10-2 


B3 


40 


0.5 


0.5 


2.5 


907 


383 


957 


459 


5.49 


10-3 


2.56 ■ 10-2 


B4 


40 


1 


0.1 


2.5 


506 


205 


438 


236 


3.04 


10-3 


1.53 ■ 10-2 


B5 


40 


1 


0.3 


2.5 


536 


216 


587 


276 


8.71 


10-3 


2.85 • 10-2 


Be 


40 


1 


0.5 


2.5 


641 


271 


676 


324 


1.55 


10-2 


7.24 • 10-2 


Br 


40 


1.5 


0.1 


2.5 


395 


179 


496 


187 


1.39 


10-2 


3 ■ 10-2 


Bg 


40 


1.5 


0.3 


2.5 


541 


185 


610 


329 


2.22 


10-2 


7.17-10-2 


B9 


40 


1.5 


0.5 


2.5 


547 


248 


602 


347 


3.2- 


10-2 


8.53-10-2 


Ci 


40 


0.5 


0.5 


5 


890 


464 


770 


349 


3.46 


10-3 


1.16 - 10-2 


C2 


40 


1 


0.5 


5 


789 


788 


772 


770 


1.35 


10-2 


8.01 - 10-3 



Table 1. Models, for different initial t^, Rq, F, and d and results for the maximum velocity, 
^maxT=timo(k™s-i), whcrc 6 is the inclination angle; and the mass-loss rate, MT=time(M0/yr). 



Mass loss within the funnel flow 

It is instructive to consider the distribution of variables within a horizontal cross-section 
at a certain height above the equatorial plane. In so doing, we interpolate the solution 
from an (r, 9) - spherical grid to a (z, w\ 100 x 100) Cartesian grid. Figure [6] shows the 
distribution of temperature and z-component of velocity, in terms of the escape velocity, 
f/esc = {'^GM^Yi/fY^'^ , at different heights for model B(^. 

A hot region extends to ^ 1 pc at ^ = 0.2, and to w ^ 2 pc at 2; = 1. The "funnel" 
can be seen in distributions of both temperature and velocity. At the X-ray heated boundary 
of this nozzle gas is being heated so that its temperature increases suddenly to ~ 10^ — 10^ 
K. This fact reveals an analogy between the torus flow with X-ray excited winds in X-ray 



binaries (IBasko et al.l 119771 : iMcCray fc HatchettI Il975l ) ; we discuss this further later in this 
section. Notice that in our case the inner surface of the torus both serves as a copious source 
of a gas and as a colhmating funnel. 



Figure [5] shows models Bq and C2 at t = 4 and Figure [8] shows density and velocity 
streamlines for models Ai and at t = 4. Notice, there is little difference between Figure 
E] (left panel, Model Ai) and Figure [8] (right panel. Model A^); the effect of smaller Rq is 
partially compensated by the fact that F is also smaller, thus reducing the effective gravity. 
If F ~ 0, then zu~ ~ 0.5 (for C = 0.2 in the equation fllll) ). However, when F = 0.5, as 
in model Bq (Figure O left panel), the effective gravity at the innermost optically thin edge 
of the torus is reduced by half. Figure [5] (right panel) shows a model with initially large 
distortion d = 5 {C = 0.1), model C2 in Table 1. 

In model Aq, a well-developed wind is observed in the vicinity of the high density torus, 
following the equal pressure contours; the maximum radial velocity is observed close to the 
axis at Wmax(^ — 3°) = 700 km s^^. As a general trend at t = 3 the maximum velocity has 
a plateau at 20° < 9 < 50°, fmax = 220 km s^^, and lower values closer to the equatorial 
plane (Figure [7]). The flow is approximately symmetric in both hemispheres. At later times, 
t = 4 and t = 5, the behavior of the model is similar to t = 3: namely, fmax(^ ^ 4°, T = 
5) = QOOkms"^, and on the plateau being fmax ~ 380km s~^. The torus is losing mass in 
all directions, although with very different speed at different inclinations. Because we are 
solving equations of ideal hydrodynamics (with only a small numerical viscosity), accretion 
through the inner boundary (at r=0) is negligible: Min(MQyr~^) < 10~^. The maximum 
mass flux per unit solid angle M^'^^(M0 yr~-^ sterrad"^) peaks at ~ 13° at t = 3, i.e. at 
much higher inclinations than Vma.x, and at 6^ ~ 55°, M^'^^ = 0.01 at t=l, M^'^'^ = 2 ■ 10~^ at 
t=3, and M^^^ = 0.02 at t = 5. The total mass-loss rate at t = 3 is M(Mq yr^^) ~ 7 ■ 10"^ 

The mass-loss rate is M(MQyr-^) ~ 2.4 ■ 10~^ at t = 4, and M(M0 yr^^) ~ 4 ■ 10"^ 
at t=5, and the change of the mass-loss rate with time is dM /dt{{MQjT~^)/jT) ~ 10^^. 
Comparing distributions of v and n we conclude, for example, that the apparent minima 
of Vp ~ Vr correlate (with a certain lag) with maxima of n and vise a versa, reflecting 
conservation of mass flux. 

As in model Aq, the model Bq funnel wind carries mass flux which doesn't change much 
during the evolution. The maximum velocity is as high as ~ 1000 kms~^ near the axis, and 
typically 200 < v^^ax ^ 600 km s^^ at 15° ^ < 50°. The bulk of the gas, which potentially 
may produce warm absorber features, moves with comparable speed. However the largest 
observed velocity in Bq model is fmax(^ — 3°) = 1200 km s~^, at t=5. The mass-loss rate is 
M{Mq yr~i) ~ 3.4 ■ 10~\ at t=3, and M(Mq yr^^) ~ 7 • 10"^, at t=5. 
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Spectra 

Computing absorption spectra is a key test for the warm absorber flow model. Several 
sample spectra are shown here, although the detailed discussion of methods and results of 
calculations of such spectra is postponed to a later paper. 

Figure [9] shows the model Aq spectrum observed at different inclinations. This figure 
shows the warm absorber spectrum at t = 3 and at t = 4. At t = 4 a rich X-ray line 
absorption spectrum exists in the range 43° ^ < 52°, and in the range 47° ^ < 55° at 
later times, t = 5. 



At t = 3 the Bq model predicts a rich spectrum for 42° ^ 6 ^ 47°. At later times 
a similar spectrum appears at lower inclinations. Figure [10] shows the model Bq spectrum 
observed at different inclinations at t = 4. At t = 5 the spectrum exists between 45° < 
6 < 50°. Notice that the region of the funnel wind in this model is bounded by the area 
unshadowed by the torus: 0° < < 40°. At > 30° column density becomes A^23 — 0.45 and 
the ionization parameter is ^ < 20. At higher inclinations the X-ray flux in the 1 < E < 2 
keV range becomes severely absorbed. 

Figure [TT] shows the evolution of the observed properties of the warm absorber flow with 
time (in the same time units) for model ^3. It can be seen that warm absorber spectra are 
changing slowly on a timescale At ~ 1. This is typical for most of our models and shows 
the range of times over which our solution can be considered as a representation of a steady 
state warm absorber flow. 

A quantitative analysis of our synthetic spectra and comparison with observations will 
be done in a later paper. This is due in part to the need for full 3-dimensional treatment of 
the transfer and scattering of line photons, which we do not present here. Rather, the spectra 
in figures [9], [10] and [11] are calculated assuming pure absorption. We can calculate crudely 
some of the properties of individual lines, and show that these are generally consistent 
with observations. A convenient way to do this is to discuss the profile of what is likely 
to be the strongest line in any synthetic spectrum, the La line of OVIII. In model Aq at 
t = 3 the full width at half-maximum (FWHM) of this line is ~ 200km s^^ . Closer to 
the BH, the maximum observed velocity is greater, i.e. models ^3 and B^ give FWHM 
~ 400km s~^ at ~ 43°, and 40°, respectively. The centroid energy of the line is at a 
blueshifted velocity (50 — 200)kms~^ with respect to line center. These velocities are less 
than th ose observed frorn, eg., NGC 3783, but are comparable to those observed from other 



objects (IMcKernan et al.ll2007l ). Such comparisons should also include the effects of scattered 
emission, which may skew the line centroid and red edge, and which we have not considered 
here. 
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Analytical estimates of the mass-loss rate 



The mass-loss rate found from numerical calculations is in approximate agreement with 
theoretical expectations. The value of the mass-loss rate, M, can be estimated by integrating 
the average mass flux (j) over the surface area of the torus exposed to X-ray radiation, 
S ~ 27r^ - Rn /A, where A = Rq/ H ~ 1. (7) may be est i mated using the same arguments as 



those of iBasko et al.l (119771 ) and iMcCray fc HatchettI (119751 ) . Namely, heating from a BH 
creates a narrow transition layer, a "skin" on the surface of the torus. There, temperature 
rises almost discontinuously from inner "cold" (T ~ 10"^, T < T^n), to outer "hot" (T > Tvir) 
value. This transition can be seen in Figure O 

Matching momentum, p + pv^ and mass flux, j = pv below and above this discontinuity, 
we obtain a well known relation: = {Ph — Po)/(Po^ ~ Ph^)^ where subscripts and h 
refer to values below and above the discontinuity. Being heated, the gas expands and its 
specific volume, V = 1/p increases. Above the discontinuity the flow is assumed to be 
isothermal so that P ~ ^/V. In the P-V plane, the transition between points Pq, Vq and 
Ph, Vh goes through the straight line with an inclination, (P/, — Po)/{Vh — Vq) > {dP/dT)T, 
and it follows < —{dP/dV)T = PhC^shi where Cs = {TZT/pY^^ is the velocity of sound. 
Since j = PhVh, it follows - Ufa < Cs^h. and the flow ini mediately above the discontinuity is 
subsonic (IBasko et al.lll977l : iMcCray fc Hatchettlll975l ). From the momentum conservation, 
Pm = Ph + Ph^l, and the mass flux associated with such heating, can be estimated as 

- . , , — TToT) where A^^ is the Mach number above the discontinuity, Pm is the 

VhMh{l + Mfj 

pressure below the discontinuity, and for simplicity we assumed Vh — Cs^h- 

McCray & Hatchett (1975) have calculated the state of the gas in the optically thin 
layer of a stellar atmosphere heated by X-rays. From their results it follows that the relation 
between Pm and Fx can be cast in the form: Pm = 10~^^a-i2- Pc, where a 12 ^ 1, reflecting 
the shape and effective temperature of the incident spectrum ( IBasko et al.lll977l ). Although 
it is essential (in order to obtain stationary transonic flow, correctly matching boundary 
conditions at infinity) that the flow above the discontinuity is subsonic, we assume that the 
sonic surface is located not far from the discontinuity, estimating V2 = Cg h, M-h = 1- Next, 



Po'- 
(j) 



we write: 



[I + A), where A is the effective X-ray albedo of the X-ray heated skin 



and we take A = 0.4 (which we simplistically assumed to be optically thin), and assume 
p = 0.5. Calculating M =< j > S, we finally obtain: 



M(M0/yr) ^ 0.16- 



(12) 



where Th^e is the temperature above the discontinuity in units of lO^K. Inserting relevant 
parameters, such as F = 0.5, /x = 0.5,R = 1 pc, Mq = 1, A = 1, and adopting the value 
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of Th taken from our numerical model Aq, Th — 10^ K, we estimate the mass-loss rate: 
M(M0 yr"^) ^ 4 ■ 10~^. Comparing results from this approximate formulae with those 
summarized in the Table 1 we conclude that they are in good accord. Given the torus mass 
in the Aq model M^or = 9.3 ■ 10^ Mq, we conclude that it may sustain such mass loss for 
~ 1 ■ lO^yr. The upper limit may be inferred from Table 1, and is found to be ~ lO^yr. 



Adiabatic loses 

The characteristic time scale at which the energy is deposited to the flow via Compton 
processes, can be cast in the form: 

T 



tx(s) ~ 9.4 ■ 10 



r^ 

10 pc 



(13) 



rf.T-i' 

where T = T /T^ and = 2.9 ■ lO^K. This should be compared with the dynamical time, 
tdyn of the flow: 



yn\ 



4.3 ■ 10^° r 



pc 



T. 



(14) 

When tdyn ^ the outflowing gas departs from thermal equilibrium and one must account 
for adiabatic losses, Aad, when calculating the temperature of the gas. Notice that the 
properties of the two-phase (o r multiple-phase) gas are conventionally described by the S- 
curve on the T — S diagram (IKrolik et al.lll98ll ). where S = Fx/{nkTc) is the other form 
of the ionization parameter. That is on the T — S plot those places where dT/dE > 
are stable to isobaric perturbations. Places where dT/dE < are unstable. Includin g Aad, 
may signiflcantly lower the temperature of the hot phase (jChelouche fc Netzeiil2005l ). This 
temperature can be estimated by equating the Compton heating rate, Aic 



4A;F.^Txto 



rripC^ 



the adiabatic losses rate, Aad ~ ~pc^. The flow near the funnel walls is less divergent than 

r 

it would be in the case of a spherically-symmetric wind, in which case the latter expression 



is a factor of 2 larger. Assuming that above the discontinuity v 



we obtain: 



ThiK) ~ 5.7 ■ 10*^ 



f.ry/^ 



'-pc 



(15) 



which gives T/j ~ 2 ■ lO^K, for parameters adopted in this paper. This value is in good 
agreement with the value of Th, which is found from T{vu, z) distributions shown in Figure 
El Three major regions within the funnel flow may be emphasized: i) a "discontinuity" 
where temperature is rising from the inner "torus" value to ~ lO^K; ii) a "plateau" where 
T ^ Tfi and thermodynamic characteristics of the flow result from the interplay between Aic 
and Aad; hi) region of hot, over ionized flow where T T^. 
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Returning current 

From Figures ([5]) and ([8]) we see that there exists a region, behind the dense torus, 
where outflow is switched to inflow. This gas rejoins the torus in the shadowed region. For 
example, taking the model Aq, and integrating the mass flux over the region where Vp < 0, 
we obtain Mjn = 4 ■ 10^^ Mq yr^-*^ at t=3. That makes ~ 6% of the total accretion rate 
MaccT, required to maintain O.SLedd luminosity of the BH, given the efficiency of accretion, 
7] = 0.06. At the same time, much more mass, ~ 2 ■ 10~^MQyr~^, is lost within the funnel 
{6 < 50°) in the X-ray excited wind. Matter that is removed from the funnel is replaced 
by gas from the torus interior. Thus, a weak large scale convection flow is observed in the 
simulations. This effect is most clearly seen in models with large d, such as model C2, shown 
in Figure [5] and is due to a strong drop of Vp as the outflowing gas is passing the shock wave 
front behind the torus (c.f. Figure [3]) and being unable to escape from the potential well. 



Radiation force 

The dependence of the radiation pressure on the ionization parameter, ^, is determined 
by equations ([H]) and fllOl) . In the region of the fast flow the wind is too overionized for 
the radiation force to be important. This resembles Low Mass X-ray Binaries (LMXB) 
case in which the radiation pres sure is also found to be insufficient to drive a significant 



outfiow (iProga &: Kallmanll2002l ). The ionization parameter drops below ~ 100 at ^ > 40° 
- the value determined by the torus aspect ratio, A. Thus the radiation pressure may be 
of importance at higher 6 and at these inclinations its relative strength is determined by 
the attenuation of the X-ray and UV fiuxes. For model Aq we have r ~ 6 at ~ 90° and 
r ~ 1 at 6* ~ 60°; i.e. the torus becomes Compton thin at < 60°. The radiation pressure 
exhibits complicated behavior with varying 6, having multiple maxima and minima. The 
force multiplier, M(t,^{9)) peaks at 44° at r ~ 2, where S'rad/s'grav ~ 5. Generally, two 
maxima of (/rad are observed at a given 6 along a radial line. The second peak becomes smaller 
at higher inclinations, i.e. in models Ai the radiation pressure is determined mainly by the 
properties of the X-ray heating (i.e. ^{0)) rather than by the attenuation of the UV fiux. At 
higher 6 smaller maxima occur at smaller r; the inner skin of the torus exerts considerable 
radiation pressure, although at large 9 it is opposed by the back pressure of the torus interior. 
We calculated a model which has the same parameter values as model Aq but with (^rad = 0. 
At t = 3 this model gives = 10°) = 564 km s"^ and v'^^''{e = 45°) = 194 km s^^ 

Comparing with Table 1 values we see that for the range of angles where warm absorber 
fiow is observed the radiation pressure doesn't play a major role in the fiow acceleration. 
In models Bi the attenuation is much stronger than in models A^ and consequently the 
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secondary maxima of (^rad which were observed in models Ai are suppressed by the e~'^ 
attenuation. The radiation pressure is important only on the skin of the torus but almost 
everywhere points in the wrong direction, opposing the back pressure of the torus interior. 
Only at 6' ~ 45° it points in the direction tangential with the torus surface, but as r ~ 3pc 
the density drops and ^ rises so that M(t,^) becomes small. 



Dependence on T, Rcfl,Tj^, and d 
If the interior of the torus is optically thick to X-rays then the torus loses mass mostly 



from the surface, much as in the 'self-excited wind' scenario for X-ray binaries (IBasko et al. 



19771 ). As shown above, in such a case the torus throat serves as a funnel and the gas is 



injected to the flow from the funnel walls. 

Notice that the location of the narrowest part of this funnel determines the characteristic 
terminal speed of the wind. In order to explore this, we have made a set of runs similar to 
models A3 and Aq, but with reduced r^. For model ^3 which has = 2 and Rq = 0.5, we 
find that for 6 = 10° and t = 4.5, the maximum velocity v^'^^ equals 738 km s~^. For model Aq 
which has = 2 and Ro = 1, we find v^'^^ = 432 kms~^ for the same 6. If in the latter model 
we make optical depth smaller, = 1 we obtain: v^^^{T = 4.5, 6 = 10°) = 400 km s^^. 
This shows, in accord with our expectations of the mass flux conservation, that the torus is 
losing mass from deeper inside. As shown in Table 1, reducing Rq has the effect of increasing 
the maximum velocity. An increase of F has the same effect. However, this maximum 
velocity may be observed at a different inclination. Increasing the distortion parameter, d 
has an effect of some increase of the maximum velocity, redistributing v^'^^{6) to higher 
inclinations. From numerical solution we notice that the torus aspect ratio, A = r/if ~ 1 
does not strongly influence the evolution. That is because it is the most inner part of the 
throat which determines the dynamics of an evaporative flow. This inner throat is located 
at high 6 so that it remains optically thick most of the time. Numerical experiments conflrm 
that the geometry of this inner throat remains approximately unchanged in time. 

Figure [2] shows that the geometry of the innermost part of the torus, i.e. the densest 
part (roughly located between 0.5 and 2 pc) shrinks considerably in the vertical direction 
during the process of the evolution. This is the result of the joint action of the radiation 
pressure and the back pressure of the hot evaporative flow. This is particularly interesting 
as it resembles the geometric ally thick outskirts of AGN ac c retion disks which are know n 



to be unstable to self-gravity (IKolykhalov fc Syunyaevlll980l : IShlosman &: Begelmanlll989l ). 



The physics of such systems is comphcated, and is subject to various possible competing 
effects. The self-gravitating instability may operate also in the torus body, perhaps leading 
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to a d ynamical system of molecular- dusty self-gravitating clouds (as in iKrolik fc Begelman 
( Il988l )). If this is the case, the optical depth of the torus, is crucial as in the optically 
thin case t he torus will effectively cool and collapse to a thin disk with subsequent star 
formation (IToomrd Il964j ). In the other extreme (r^ >> 1) the released energy can go 
to increase the v elocity dispersion of the clouds, effectively supporting the torus thickness 
(jPaczynskil 1 19781 ). Strong IR radiation pressure exerted on these clouds, which can come 



from internal reprocessing of X-rays, can produce significant vertical force (IThompson et al. 



20051 : iHonig fc BeckertI 120071 ) . and may suppress the sel f-gravity inst ability and at the same 
time provide pressure support against vertical collapse (IKrolikl 120071 ). Vertical support, and 
partial suppression of gravitational collapse, may also b e provided by radiation pressure 
from star formation within the torus or the obscuring flow (IWada &: Normanll2002l ). Further 
heating and loosing mass induces a torus to expand and change of shape. 

We have calculated models Bq and Aq with 100x100 resolution further in time to learn 
the late time evolution. At T = 17 in the Bq model the torus has two extended lobes in 
both hemispheres with an opening angle 45°. They have a certain degree of asymmetry with 
respect to the equatorial plane. The shape of the obscuring structure no longer resembles the 
initial torus; the column densities are in the range of N23 = 10 at 30°, and A'"23 = 10^ at 30°. 
The radial velocities in this structure are in the range of 200 — 400 km s"^. Our model does 
not allow for the replenishing of the torus; obviously the torus will evaporate completely if 
given enough time. Thus, in the the model Aq the torus evaporates completely by the time 
T = 15. These results imply that in order to get a quasi- stationary warm absorber flow the 
replenishing time should be of the order of the mass-loss time. The whole torus configuration 
may be unstable in a secular sense; the instability is driven by the long characteristic time 
of the global torus heating/cooling (due to expansion, winds, radiation loses), advection of 
heat in the torus body by internal flows etc. For example, the mass-loss rate, M ~ E; the 
surface area E increases during the torus expansion. If after some time of extensive heating, 
the torus separates into several parts, further mass loss will increase due to the larger total 
surface area of the fragments. 



6. Conclusions 

We have studied X-ray excited winds from the putative gas-dusty torus in AGN. We ap- 
proach this problem using numerical methods combining detailed hydrodynamical modeling 
with calculation of the warm absorber spectra. Hydrodynamical calculations include two- 
dimensional, axially-symmetric rotating flow, driven primarily by X-ray heating. Compton, 
bremsstrahlung, and photoionization heating/cooling processes were taken into account as 
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well as the radiation pressure force, which was calculated in the Sobolev approximation. A 
code combining XSTAR for photoionization calculations with the Sobolev radiation transfer 
has been developed for the calculation of the spectra. 

We find that a rotationally supported torus heated by radiation from the inner accretion 
disk and black hole can indeed be a source of the material we observe in the warm absorber 
flow. Wc find that the inner throat of the torus is not only important as a source of the gas 
but also because it creates a funnel for the outflowing wind. This leads generally to larger 
velocities within the funnel, and different velocity distribution within the warm absorber flow 
from those derived from models based on spherically-symmetric winds. The wind mass-loss 
rate within the funnel is not very sensitive to the details of the initial torus distribution and 
approaches ~ 0.02 — O.O9M0yr~^. Strong X-rays heat the gas within the funnel, producing 
a fast, ~ 1000 km ionized flow near the axis, and slower, < 500kms~^, flow closer to 
the funnel walls. This is where optical depth effects become important and a warm absorber 
spectrum is produced. Using methods developed in studies of X-ray binaries we were able 
to estimate the mass-loss rate from such funnel flow, flnding it to be in a good agreement 
with our numerical solution. 

The funnel flow is found to be promising with respect to obtaining high velocity warm 
absorber flows. What is beyond the scope of our models is the possibility of having mul- 
tiple phases in such high velocity flow, on spatial scales smaller than our grid resolution. 
Our treatment of the gas thermal properties will produce two-phase behavior at our grid 
resolution; we do not find this behavior, owing to the fact that the cooling timescales are 
generally too long. The answer to the question of whether there can be high velocity 'bullets' 
or 'embedded clouds ' on length scales smaller than the resolution of the grid is related to 
the problem of the origin of broad and narrow UV/optical line emitting clouds, and requires 
different computational methods from those employed here. 

Our models which have initial Compton depths ^ 1, aspect ratio Rq/H ~ 1, and 
located at 0.5 ^ r < 1.5 pc predict warm absorber spectra, thus confirming the main 
conclusion made in paper 1. The existence of such spectra depends on the fact that the 
fiow is intrinsically two-dimensional, meaning that both the dynamics of the funnel fiow is 
different from ID models and optical depth effects are important as they strongly depend on 
inclination. The latter point requires that we include the entire torus in the computational 
domain rather than considering it as a boundary condition. The distribution of the ionization 
parameter, ^ depends strongly on 6, further confining the range of angles where conditions 
are right for the warm absorber fiow to be observed. In most of our models warm-absorber- 
like spectra are produced in a 10° range, at ^ ~ 40 ± 5°. This range is set both by the 
initial aspect ratio of the torus, which we take to be ~ 1, and by the thickness of the X-ray 
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heated 'skin' of the torus. More optically thin models produce warm absorber-like spectra for 
^ ^ 40 ± 10° , as they potentially provide more partially optically thin gas for evaporation. 

The bulk of the gas in this scenario has a terminal velocity of the order of the escape 
velocity at the inner torus edge. Because of the funnel mechanism part of the gas is re- 
distributed to lower inclinations and acquires a higher terminal speed, ~ 1000 km s~^. In a 
real AGN environment such flow may contain clumps and irregularities and even dust, which 
are not captured in our studies because of the intrinsic limitations our methods. Accounting 
for the multiple phases of a gas (on a subcellular level) may reveal this in more detail and 
may also broaden the range of angles where the warm absorbers appear. 

The part of the flow that is shielded by the optically thick part of the torus body can 
also flow out as part of a torus global expansion. Thus it strongly depends on the deposition 
of energy directly to its interior. This problem is related to o ne of the infra red support of the 



AGN torus vertical structure against gravitational collapse ( iKrolikI 120071 ) and also requires 
additional investigation. 
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LOG DENSITY 




x (parsecs) 

Fig. 2. — Model Bq; Color-intensity plots of the logarithm of the dimensionless density 
a function of time which is given in orbital periods. In the northern hemisphere this 
superimposed with velocity vectors. Axes: distance in parsecs. 
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Fig. 3. — Model Bq, Color- intensity plots of the logarithm of the dimensionless pressure at 
t = 3; at northern hemisphere superimposed with velocity vectors. The location of the sonic 
surface is marked in red. Axes: distance in parsecs. 
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Fig. 4. — Model Bq; distributions of Vr, and the poloidal velocity Vp (upper left), ^ (upper 
right), n (lower left), and the rate of growth of number density with radius (lower right) 
at time, t = 3. Curves are marked by a inclination angle 6. Horizontal axis: distance in 
parsecs. 
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Fig. 5. — The effect of distortion parameter d; Velocity streamlines superimposed on contours 
of the number density. Model Bq (left panel), model C2 (right panel) at time, t=4. Axes: 
distance in parsecs. 
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Fig. 6. — Model Bq at time, t = 3. Horizontal 'slices' of temperature and z-component of 
velocity are marked at each curve by the value z in parsecs. Vertical axis; left: z-component 
of the velocity in terms of local escape velocity; right: temperature. Horizontal axis: distance 
in parsecs. Curves; dashed: velocity; solid: temperature. 



- 30 - 




Fig. 7. — Model Aq; As in Figure H] here are shown the distributions of Vr, Vp (upper left), 
C, (upper right), n (lower left), and the rate of growth of number density with radius (lower 
right) at t = 3. Curves are market atop by an inclination angle 6. Horizontal axes: distance 
in parsecs. 
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Fig. 8. — The effect of F; velocity streamlines superimposed on contours of the number 
density. Model Ai (left panel), model A^ (right panel) at time, t=4. Axes: distance in 
parsecs. 
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Fig. 9. — Model Aq. spectrum, observed at time, t=3. 




Fig. 10. — Model Bq: spectrum, observed at time, t=4. 
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Fig. 11. — Model ^43: X-ray spectra, observed at ^ = 45°, as a function of time. 



